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The recently fabricated two-dimensional magnetic materials CugX2(cpa)6 • sEbO (cpa=2- 
carboxypentonic acid; X=F,Cl,Br) have copper sites which form a triangular kagome lattice (TKL), 
formed by introducing small triangles ( "a-trimers" ) inside of each kagome triangle ( "b-trimer" ) . We 
show that in the limit where spins residing on b-trimers have Ising character, quantum fluctuations 
of XXZ spins residing on the a-trimers can be exactly accounted for in the absence of applied field. 
This is accomplished through a mapping to the kagome Ising model, for which exact analytic so- 
lutions exist. We derive the complete finite temperature phase diagram for this XXZ-Ising model, 
including the residual zero temperature entropies of the seven ground state phases. Whereas the 
disordered (spin liquid) ground state of the pure Ising TKL model has macroscopic residual entropy 
ln72 = 4.2767... per unit cell, the introduction of transverse (quantum) couplings between neigh- 
boring a-spins reduces this entropy to 2.5258... per unit cell. In the presence of applied magnetic 
field, we map the TKL XXZ-Ising model to the kagome Ising model with three-spin interactions, 
and derive the ground state phase diagram. A small (or even infinitesimal) field leads to a new 
phase that corresponds to a non-intersecting loop gas on the kagome lattice, with entropy 1.4053... 
per unit cell and a mean magnetization for the b-spins of 0.12(1) per site. In addition, we find that 
for moderate applied field, there is a critical spin liquid phase which maps to close-packed dimers 
on the honeycomb lattice, which survives even when the a-spins are in the Heisenberg limit. 



PACS numbers: 75.30.Ds, 75.10.Hk, 71.27.+a 



I. INTRODUCTION 

Geometrically frustrated spin systems hold promise for 
finding new phases of matter, such as classical and quan- 
tum spin liquid ground states. Of considerable interest 
has been the discovery of a stable phase with deconfined 
spinous in a model of quantum dimers on the (geomet- 
rically frustrated) triangular lattice, i.e. a spin liquid. 1 
Beyond the interest in fundamental theoretical issues of 
these modelsji^ physical realizations of these systems 
may have technological applications in achieving lower 
temperatures through adiabatic demagnetization. Since 
such techniques require a material which can remain in a 
disordered, paramagnetic state to very low temperatures 
(rather than undergoing a phase transition to an ordered 
state), this makes geometrically frustrated spin systems 
attractive for such applications. 

Recently, a new class of two-dimensional magnetic ma- 
terials CugX2(cpa)6-xH20 (cpa=2-carboxypentonic acid, 
a derivative of ascorbic acid; X=F,Cl,Br) — was fab- 
ricated in which Cu spins reside on a triangular kagome 
lattice (TKL), formed by inserting an extra set of trian- 
gles (a-trimers) inside of the kagome triangles (6-trimers). 
(See Fig. Q]) In a recent paper—, we analyzed the ther- 
modynamic behavior of Ising spins on this lattice using 
exact analytic methods as well as Monte Carlo simula- 
tions in finite held. In this paper, we extend our analysis 
to include quantum fluctuations of the spins on a-trimers, 
i.e. we study an XXZ-Ising model on the TKL. The Cu 
spins in CugX 2 (cpa)6 ■ xH 2 likely have isotropic Heisen- 
berg interactions, for which exact solutions are currently 



inaccessible on a frustrated lattice. However, in the limit 
where spins on 6-trimers have Ising character, it is still 
possible to take into account the quantum fluctuations 
on a-trimers exactly through a mapping to the kagome 
Ising model. In the presence of applied field, the model 
maps to the kagome Ising model with three-spin interac- 
tions. We present exact results for the phase diagram at 
all temperatures without applied field, and at zero tem- 
perature in the presence of applied field. 

In the absence of applied field, the zero temperature 
phase diagram is richer than the case where both a and 
b spins are in the Ising limit. The ordered phase sur- 
vives quantum fluctuations, but for large enough trans- 
verse coupling, there is a first order transition to a state 
with lower total spin on the a-trimers. For the disor- 
dered phase, the quantum fluctuations of the a-spins 
partially lift the ground state degeneracy. In the pres- 
ence of applied field, the phase diagram is even more 
rich. The "honeycomb dimer" phase, which is present in 
applied field for the purely Ising version of the model^ 
survives the introduction of quantum fluctuations. How- 
ever, rather than arising for infinitcsimally applied field, 
a finite field is now required. Infinitesimal applied field 
in the presence of quantum fluctuations of the a spins 
results in a new phase, which we have mapped to a non- 
intersecting loop gas on the kagome lattice. 

This paper is organized as follows. In Section [II] we in- 
troduce the XXZ-Ising model on the TKL. In Section Hill 
we present an exact mapping to the kagome Ising model, 
and we derive the finite temperature phase diagram in 
zero field. In Section [TVl we study the finite field case, 
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FIG. 1: (Color online) The triangular kagome lattice (TKL) 
with the nearest-neighbor couplings J aa and J z b . Solid circles 
are "a" sites and open circles are "b" sites. The shaded region 
represents the unit cell. 
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FIG. 2: (Color online) A hexamer consisting of three a-spins 
(filled circles) and three 6-spins (open circles). The coupling 
between a-spins is of the XXZ type, whereas the a-b coupling 
is of Ising type. The mixed Heisenberg-Ising model on the 
TKL can be fruitfully analysed by breaking it up into hexamer 
units. 



using an exact mapping to the kagome Ising model in a 
field with a three-spin coupling term. In Section [V] we 
present our discussion and conclusions. 



II. XXZ-ISING MODEL ON THE TKL 

The TKL can be obtained by inserting an extra set 
of triangles inside of the triangles of the kagome lattice. 
Alternatively, it can be derived from the triangular lat- 
tice by periodically deleting seven out of every sixteen 
lattice sites. This structure has two different spin sub- 
lattices, a and 6, which correspond to small trimers and 
large trimers, respectively. Each spin has four nearest 
neighbors. The unit cell contains a total of 9 spins (6 on 
the a sublattice, and 3 on the b sublattice). The space 
group of the TKL is the same as that of the hexagonal 
lattice, p6m, in Hermann-Mauguin notation. The shaded 
region in Fig. Q] encompasses one unit cell of the TKL. 

The Cu spins in real materials have S = 1/2, and quan- 
tum effects cannot be neglected a priori. The Ising limit 
of this model has been previously considered by Zheng 
and Sun£, as well as by three of us£. The ground state 
phase diagram^ as well as many experimentally testable 
thermodynamic quantities^ have been calculated. In this 
paper, we include the quantum fluctuations of the a spins 
(i.e. those on a-trimers), while treating the b spins as 
classical Ising spins. In this limit, the a spins can be in- 
tegrated out exactly, leaving an Ising model of the b spins 
on a kagome lattice, which we solve exactly. 

We consider a model in which the exchange coupling 
between neighboring a-spins is of the XXZ type and the 
coupling between neighboring a- and 6-spins is Ising-likc. 
Whereas this treats the a-spins as fully quantum mechan- 



ical, the 6-spins are classical. The Hamiltonian is 

h = - £ [J z aa sfsf + j^{sfsf + sYsJ)] 

J2 Ja b sfsf -h^s: a) 

(iea.jeb) i 

where Sf' Y ' Z are the S = 1/2 spin operators at site i, 
angle brackets indicate summations over nearest neigh- 
bors, and h is an external magnetic field. With this sign 
convention, positive coupling J > corresponds to fer- 
romagnetic interactions, and negative coupling J < is 
antifcrromagnctic. This model contains the following en- 
ergy scales as parameters: J^ a , J* a , J^ b , T, and h. We 
will take | J^ b | as the unit of energy. 

When J*,, = 0, this model reduces to the Ising model 
on the TKL, which was studied exhaustively in Ref. 0. 
(See also Ref. d.) The J^ a terms permit transverse quan- 
tum fluctuations within each a-trimer. However, the 
quantum fluctuations are confined to reside within the 
a-trimers, as can be seen by the fact that the total S z of 
each a-trimer is a conserved quantity. 

In spite of the frustrated nature of the model, it turns 
out that it can be exactly solved at all temperatures for 
h = and at zero temperature for finite h. The Hamilto- 
nian can be written as a sum over hexamers, (see Fig. [5]) 
where each hexamer consists of an a-trimer and its en- 
closing 6-trimcr, H = J2 n Hn- The Hamiltonian for each 
hexamer is given by 

— Jab^ai^bj ~~ ^ &ai — ^ ^bj ■ (2) 

i j^Li i j 

Note the factor of in front of 5^ , which is required 
because each 6-spin is shared by two hexamers. 
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Different H n commute with each other. Furthermore, 
each a-spin appears in one H n but no others. Therefore, 
we can perform the trace over all of the a-spins to give an 
effective Hamiltonian involving only o-spins as follows: 



Z = Tiexp(-/3H) 
= > Tr exp 

{S,} 



-f3H(S b ,S a ) 



(3) 
(4) 

(5) 



The contribution to the partition function coming from 
the trace over the a-spins within each hexamer depends 
on the values of the surrounding 6-spins, i.e.. 



The trace can be evaluated by diagonalizing the 
hexamer Hamiltonians for all eight configurations of 
the enclosing o-spins, in order to obtain the en- 
ergy eigenvalues Ej{S bl , S b2 , S b3 , h), j = 1,2,3, 8 
(see Fig. [3j), and subsequently calculating Z({S b }) = 
Y]j exp [—PEj({S b })]. The energy diagonalization is par- 
ticularly simple since the total S z of each a-trimer is 
a good quantum number. In the case where the sur- 
rounding 6-spins are in the (TTT) or (III) configura- 
tions, total >S 2 is also conserved. Due to the local C 3 
symmetry of each hexamer, Ej(\H,h) = £j(T|T,/i) = 
Ejdl^h). In addition, the energy eigenvalues re- 
spect time-reversal symmetry, so that Ej(Hl,h) = 
Ej{m,-h), and E j (Ul,h)=E j (M,-h). 



Z(S bl ,S b2 ,S b3 ) = Tr exp 

Sal ,S a 2-,S a 3 



-f3H n (S b ,S a ) (6) 





4^-CTTi, ft) 


— ^Jaa + Jaa ~ U ah — 5ft 






2J a a T ^aa ^Jab 3ft 




2Jaa + Jaa ~ 2J a() — 5ft 






— 3Jaa ~~ 2Jab ~ 7/l 




2Jaa T -Jaa ~ ^J ab ~ 5ft 






— 3J aa + 5ft 




-3J* - 6 - 9ft 






2Jil + + 2J,f b + ft 




-3J* + 6Jf 6 + 3ft 




+ J z 




+ ft 


-iJaa + Jaa + 2J& - ft 




~ ''aa 


+ V-4^a X a+9Ja Jf a 2 +4J^ 


+ ft 


^Jaa + Jaa H~ ^^ab ~ ^ 


-J x 


1 aa 


-V / 4^^ + 9Jf 2 +4J V- 


- 3ft 


2</ a a H~ 'Jaa H" ^Jab ^ 


-J x 

aa 


+ J Z 

1 aa 


+ V4^^+9^ 2 +4J^, 2 - 


-3ft 



TABLE I: The eight energy levels of a hexamer arising from the a-spin degrees of freedom, for two configurations of the enclosing 
6-spins (HI and TTT) - 



Hence we find 

_|_ e 4( b ' l + 4J aa- J aa+ ZJ a b)P _|_ g 4 4 2 _|_ g T 1 4 1 2 (7) 

Z(TTI, /i) = e-i( 5h - 3J °°> +2J ^ + ei {3h - 2J ^ J ^ +2J -" )p + e -^ h+2J ^+ J L+ 2J a b )0 + e i(^+3Jf a +2J^ 

e J (-h+ J* - Jaa+V 9J ^ 2 - iJ ab J a a +^Sb) , g~ 1 J fa + J fa + V 9J ^ 2 -4 Jf b J£ +4 Jf b 2 ) /3 

+ e ^(3/ i +J; Y a -J a z a + V9J„ A ; 2 +4J a z 6 ./ ,f £l +4./f b 2 )/3 + e - i (-3h- J* + J Z + ^9 J a x Q 2 +4 Jf„ j£ +4 Jf^) , g . 

zau,ft) = z(m,-ft) (9) 

z(TU,*0 = -Z(m,->0- (io) 



III. ZERO FIELD kagome lattice exactly, so that up to a temperature- 

dependent additive constant, we have 

A. Exact mapping to the kagome Ising model 

When h = 0, the trace over a-spins maps the XXZ- H c s({<r b }) = —J bb auCbj, (11) 

Ising model to the ferromagnetic Ising model on the (ij) 
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where Jtb is the effective 6-spin coupling and au = ±1 
for consistency with the Ising model literature. 

The simple form of H c s (two-spin nearest neighbor in- 
teractions only) can be understood easily. Since the hex- 
amer Hamitonians commute between different hexamers, 
the trace over a-spins in a hexamer couples only the three 
6-spins in that hexamer. When h = 0, H c g must have 
global up-down symmetry and so it cannot contain odd 
powers of a : Taking into account that a 2 = 1 then im- 
plies the quadratic nearest neighbor H g above. To find 
Jbb, we evaluate the two different actions Z(TTt) an d 
Z(TTI) (others are related by permutation and up down 
symmetry of the 6-spins) and match them to 



Z(&bi, a b2i Vbz) = Z a exp (3Jbb(cr b2 ab3 + o- fc3 c7 h i 
whereupon we obtain 

z a = z(m) 1/4 ^UH) 3/4 , 
_i £0Tt) 

pJbb --^zaur 



C61C62) 

(12) 



(13) 
(14) 



Thus, we see that integrating out the a-spins gives rise 
to effective couplings between the fo-spins, which corre- 
sponds to the classical Ising model on the kagome lattice. 
The effective Ising coupling J b b is a complicated func- 
tion of the original couplings and J^ b , described 
by Eqs. ([T])-([I0]) and Eq. (fl"4"|) . and of the temperature. 
However, in all cases, it is non-negative and hence the 
effective 6-spin coupling is ferromagnetic. 

The ferromagnetic kagome Ising model has an exact 
solution, which has been known for some time^. There- 
fore, the XXZ-Ising model on the TKL can also be solved 
exactly. In particular, the critical temperature, free en- 
ergy, energy density, entropy, and specific heat can be 
calculated in the same manner as in Ref. 0, which we 
outline in the next section. 



B. Critical temperature, free energy, and entropy 

The kagome Ising model has a phase transition to 
a ferromagnetic ordered state at a coupling strength 



/^kagome = 3ln(3 + VT2) = 0.466566.... Therefore 
the TKL XXZ-Ising model has a phase transition when 

WLJ^Jab,® = /^kagome- Th c ferromagnetic 
state of the kagome model corresponds to a /erromagnetic 
phase in the TKL for J^ b > (since then the a-spins are 
aligned with the 6-spins), and to a /errimagnctic phase 
for J^ b < (since then the a-spins are antialigncd with 
the 6-spins). 

For convenience, we define the free energy as F = In Z . 
Based on the discussion above, we see that the free energy 
of the TKL XXZ-Ising model is given by the sum of the 
free energy of the kagome Ising ferromagnet and a term 
that arises from integrating out the a-trimers, i.e. 
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FIG. 3: (Color online) Energy levels of the hexamers when 
the enclosing b spins are in the (TTT) configuration (blue lines) 
and when they are in the (ITT) configuration (red lines). The 
energies are shown as a function of the transverse coupling 
J^, for \Jat\ = 1 and for two different values of J^ a - (The 
energy levels are independent of the sign of Jab-) The thick 
blue-red-blue lines represent triply degenerate states, which 
include two degenerate states from the (TTT) subspace of the 
6-spins, and one state from the (ITT) subspace. Panel (a) 
shows a representative phase transition from phase 7(1/2) to 
Phase IX for J aa < —\Jab\- Panel (b) shows a representative 
set of phase transitions from phase 7(1/2) to 7(3/2) to phase 
IX for J* >-|Jf 6 |. 



f( J Za> J aa, J ab, P) = MPJbb) + 2/ Q 



(15) 



where f3Jbb{ Jaa: Jaai Jab> 0) ^ s tne effective kagome cou- 
pling from Eq.Q3] The free energy per TKL unit cell is /; 
fb is the free energy per kagome unit cell; and f a = In Z a 
is the free energy contribution per a-trimer. The TKL 
unit cell corresponds to one kagome unit cell, and it con- 
tains two a-trimers, hence the factors in the above equa- 
tion. 



Other thermodynamic quantities can be obtained by 
differentiating / with respect to (3. It will be convenient 
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to define the quantities 

U j ({<rb}) = jJ2 E ^ a »V e ~ 0Ej{{ ' Th}) > 



(16) 



C j ({<r b })=0 2 U^ma^fe- 



(17) 



where {ct&} =TTT or ft J,. Then, the energy per unit cell 
is 



dp 
C/ TTT 



3U 



Til 



u 



c/ TTT 



2 4 
and the heat capacity per unit cell is 



Uka.g((3Jbb), 



du 

C TTT 



3C 



1 IJ 



C 



TTi 



ill 



+ 



\ 4 



{[3Jbb), 



(18) 
(19) 

(20) 

(21) 



kau, 



where Uk ag is the energy per unit cell of the kagome lat- 
tice Ising ferromagnet (in units of the kagome coupling) 
and Ckag is the heat capacity per unit cell of the kagome 
lattice.— i 10 ! 11 ! 12 The entropy per unit cell is s = f + flu. 
The ground-state entropy can be calculated by taking the 
limit /3 — > oo and observing that Uj({<Tb}) an d Cj({<Jb}) 
are dominated by the lowest energy levels of the hcxamcr. 
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\ J x = -2J Z - 2 IJ Z 1 

\ J aa ' LJ aa \b 



FIG. 4: (Color online) Exact T = phase diagram of the TKL 
XXZ-Ising model in the {Jaa) Jaa) plane, with \ J%b\ = 1 as 
the unit of energy. The black lines indicate phase boundaries. 
Phases 1(1/2) and 1(3/2) are described in the text; they are 
ferromagnetic if J^b > and ferrimagnetic if J^b < 0. 




C. Zero-temperature phase diagram and 
ground-state properties of each phase 

In the absence of h, the partition function is invariant 
under a change of sign of J^ b (because this can be gauged 
away by redefining the Ising 6-spins), so the key physics 
is independent of the sign of Jj, . However, because of 
the frustrated geometry of the trimer, it is not invariant 
under a change of the sign of J^ a ■ 

The behavior of the kagome Ising ferromagnet is de- 
termined by the dimensionless coupling K = pjbb- For 
T = 0, we need the value Kq which is the limiting value 
of K as P goes to infinity. In that limit, each a-trimer is 
restricted to its ground state for a given configuration of 
the surrounding 6-spins, and Eq. (Ti"4"|) becomes 

pJbb = i[hi^i + /3(£o(TU) - Sq(TTT))] (22) 

where So(TTT) an d -Eo(ttl) are the a-trimer ground state 
energies for the two relevant configurations of the 6-spins, 
and D\ and D 2 are their respective degeneracies. Pro- 
vided -Eo(Ttl) is greater than EodTt), then clearly Kq 



FIG. 5: (Color online) Finite-temperature phase diagram of 
the TKL XXZ-Ising model, with |jf 6 | = 1 as the unit of 
energy. In the dark blue region of parameter space (phase 
IX), the system is disordered even at zero temperature. The 
straight lines in Fig. [4] corresponding to phases V and VIII 
turn into critical surfaces in Fig. [5] However, the boundary 
line between 1(1/2) and 1(3/2) in Fig. [4]is hidden underneath 
the critical surface in Fig. [5] 



goes to infinity, and the 6-spins will have perfect fer- 
romagnetic order. On the other hand, if £^o(TTI) = 
^o(TTT), Ko will be finite but can still be > if D x > D 2 . 
In this case, whether the 6-spins have LRO depends on 
whether Kq exceeds the critical coupling of the kagome 
Ising model. Note that here the coupling between the 
6-spins is caused by maximizing the ground state degen- 
eracy. 

By studying the energy levels -Eo(ftl) an d i?o(TTT) 
for different combinations of parameters, we can work 
out the entire zero-temperature phase diagram in the 
{Jaai Jaa) plane, which is shown in Fig. |U Further- 
more, the values of parameters at which the energy levels 
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0.20 r 



2.5 



FIG. 6: (Color online) The dimensionless effective Ising two- 
spin coupling between adjacent 6-spins without magnetic field, 
/3Jbb, as a function of temperature T. There are four values 
of the transverse coupling J^ a = 0, —0.1, —1, —10 (bottom to 
top). Other parameters are J^ a = —3 and J^b — ±1- 



Ej (TIT) or Ej (ITT) undergo level crossings are also of im- 
portance, as it changes the value of Di and/or D 2 - We 
label the phases I, V, VII, VIII, IX for consistency with 
our previous work£. Table |TT] shows the residual entropy 
of each of these phases, calculated using the approach 
described in the previous section. 



Phase 


PJbb 


so 


so/9 




I 


CO 








(LRO) 


V 





In 72 = 4.2767 


0.47519 





VII 


Mn^ 

4 111 3 


4.4368 


0.49297 


finite 


VIII 


|ln3 


2.8918 


0.32131 


finite 


IX 


±ln2 


2.5258 


0.28064 


finite 



TABLE II: Properties of the various T = phases marked 
on the phase diagram in Fig. [4] Here, /3Jbb is the effective 
Ising coupling between 6-spins, so is the entropy of the entire 
system per TKL unit cell, so/9 is the entropy per site, and 
£ is the correlation length of the fa-spins. Comparing /3Jbb to 
the pJ c of the kagome lattice shows that only phase I has 
long-range order. 

Note that the entropies satisfy the inequalities 
Sy , Sy ii j Syiii > Six, Si and Syn > Sy, Syin- This 
agrees with the intuition that when the system is tuned 
to a phase transition line or point, the system is able 
to access states from both adjacent phases, and there- 
fore the entropy is higher than that of the surrounding 
phases. 

Fig. [5] shows the entire finite temperature phase dia- 
gram as a function of Jf /|JfJ, J*/\J* b \, and T/\J? b \. 
Apart from the 1(3/2) to 1(1/2) ground state transition, 
which corresponds only to a change in the local physics, 
it is clear that all of the phase transitions survive at fi- 
nite temperature. Not surprisingly, as temperature is 
increased, the disordered phase becomes a larger part of 
the phase diagram. 



Based on experimental data on the TKL magnets 
CugX2(cpa)6 • xH^Oj^^, it is sensible to make the fol- 
lowing assumptions: J z a and J* a are antiferromagnetic 
(JLJaa < 0), and |J*| > |J*|. At zero applied field, 
this would put these materials in the disordered phase. 



D. Physical explanation 

We now discuss the various phases and the transi- 
tions between them. Phase I has £^ (TTI) > -E-o(TH), 
as is shown in the right side of Fig. [31 and therefore 
R'o = K(T -> 0) is infinite (see Sec. H1TC]) . In this 
case the 6-spins arc perfectly fcrromagnetically ordered 
at T = 0. This ferromagnetic phase of the 6-spins is 
further subdivided into two T = phases by the line 
Jaa = Jaa + l^abl ( see Fig- E]), corresponding to differ- 
ent a-trimer configurations. In the phase 1(3/2), the a- 
trimers are in a unique ground state with total spin value 
S = 3/2 and total z-component of spin S z = 3/2. In the 
phase 1(1/2), the a-trimers are in a unique ground state 
with total spin value 5 = 3/2 and total z-componcnt 
of spin S z = 1/2. On the phase boundary, these two 
states become degenerate, and each a-trimer is two-fold 
degenerate, with an associated T = entropy of In 2. At 
finite temperature, the transition from 1(3/2) to 1(1/2) 
becomes a crossover, since the difference between phases 
1(3/2) and 1(1/2) is the local configuration of each a- 
trimer. This crossover, which is evident in the region 
Jaa = Jaa + \^ab\ m Fig- HI should be accompanied by a 
relatively sharp peak in the entropy near the T — phase 
boundary. From the figure, it can also be seen that the 
1(3/2) phase is more robust than the 1(1/2) phase, in that 
it has a higher transition temperature. This can be qual- 
itatively understood as follows. In phase I, the effective 
coupling J b at T = is determined by ^(TTI) - £b(ttt)- 
Far into phase I (large | J a a\), the a-trimer ground state is 
locked in its J z b = configuration, and _Eo(TTI)^-E'o(TTT) 
can be calculated by perturbation theory to be \J z b \ for 
1(3/2) and §|jf b | for 1(1/2). Thus, 1(3/2) has a larger 
Jbb than 1(1/2) and therefore a larger transition temper- 
ature. 

Let us consider going from phase 1(1/2) to phase IX 
across the phase boundary V. Figure [3] shows the eight 
energy levels of an a-trimer, for the two symmetry- 
distinct configurations of the 6-spins (TTT and J.tt)i with 
J z a = — 3 and J z b = dbl. The most important feature of 
the graph is how £o(TTI) com P a res with £ (TTT); along 
with their respective degeneracies. Starting from phase 
1(1/2), as J* a is decreased towards 0, the difference be- 
tween -EulTTi) and _Eo(TTT) decreases until they become 
equal at the full Ising limit (phase V), J* a = 0. The 
ground state of the a-trimer is triply degenerate in this 
limit, regardless of the configuration of the 6-spins. Be- 
cause of this, the effective coupling between 6-spins goes 
to zero, A'q = 0. This was noted in our previous paper- 
on the Ising limit of the TKL model, along with the cor- 
responding ground state entropy of In 72 per TKL unit 
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cell in the Ising limit. 

Now turn on an antifcrromagnctic transverse coupling 
J* a < to go into phase IX. The lowest energy remains 
independent of the 6-spin configurations. Notice, though, 
that the lowest eigenvalue of the (|TT) subspace is dou- 
bly degenerate, whereas the lowest eigenvalue of the (J. ft) 
subspace is nondegenerate. In the (ttt) subspace of the b 
spins, the lowest a-trimcr energy state in this region has 
5 = 1/2, S z = 1/2, and there are two independent states 
with these quantum numbers. They are degenerate be- 
cause the surrounding (ttt) 6-spins couple to the a-spins 
through S z only. When the 6-spins are in the (||T) con- 
figuration, this degeneracy is broken. This means that at 
zero temperature, the (ttt) configuration of the 6-spins 
is twice as likely to occur as the (|tt) configuration. The 
zero temperature effective coupling constant between b- 
spins is therefore Kq = \ In 2. 




FIG. 7: (Color online) The dimensionless effective Ising two- 
spin coupling between adjacent 6-spins with a finite magnetic 
field h — 0.1, /3J(,(,, as a function of temperature T, for four 
values of the transverse coupling J* a = 0, —0.1, —1, —10 (bot- 
tom to top). Other parameters are J^ a = — 3 and J^ h = ±1. 



To get a fuller picture of what is happening, we plot 
the effective coupling (3Jbb as a function of temperature 
in Fig. [6l We see that as long as quantum fluctuations 
are present (i.e. J^ a is finite), the value of f3Jbb increases 
monotonically as T is lowered, approaching a constant 
value ±ln2 = 0.173287... as T -> 0. This value is less 
than the critical coupling of the kagome Ising ferromag- 
net, /3J b c b = ±ln(3 + \/l2) = 0.466566.... Therefore phase 
IX is disordered at all temperatures. The entropy in 
Phase IX is 2.5258... per unit cell. It exists precisely 
at h = and J* a = Jaa < — 1- Note that in deviating 
from the Ising antifcrromagnctic limit (i.e. Phase V) by 
adding quantum fluctuations, the entropy decreases. 

Next we consider crossing the phase boundary from 
1(3/2) to IX through VIII. Along the line VIII, we again 
have £o(Tn) = £o(TH)- The difference from IX, how- 
ever, is that when the 6-spins are in the (ttt) configura- 
tion, the 5 = 3/2, S z = 3/2 state of the a-trimer is also 
degenerate with the two 5 = 1/2, S z = 1/2 states. This 
leads to a larger value of Kq = j In 3, and hence a longer 
T = correlation length, than phase IX. Nevertheless, 
this value of Kq is still below that of the critical value of 
the kagome Ising model, and the ground state remains 
disordered. 

If we are in phase IX close to the line phase VIII, 
the a-trimer ground state when the surrounding 6-spins 
are all up is doubly degenerate, with a small gap to the 
5 = 3/2 state. If the temperature is now increased, then 
this low-lying excited state will have non-zero Boltzmann 
weight, and hence the value of K will first increase be- 
fore decreasing again at higher temperature due to the 
effects of thermal disordering. Thus, the change oiK and 
hence the o-spin correlation length with temperature will 
be non-monotonic, first increasing with temperature be- 
fore decreasing. This is a manifestation of the order-by- 
disordcr mechanism commonly seen in frustrated systems 
with residual ground state entropy. 



IV. FINITE MAGNETIC FIELD 

A. Exact mapping to the kagome Ising model with 
three-spin interactions 

We now consider the XXZ-Ising model on the TKL 
in the presence of finite magnetic field. This model is 
defined in Eq. ([T]). Note that we consider applied field 
parallel to the axis of the Ising spins on the b sublattice, 
so that the b spins remain classical in their behavior. In 
the presence of finite field h, spins on the a-sublatticc 
can still be integrated out, yielding an effective model 
in terms of the b spins only, which reside on a kagome 
lattice. However, since the original model for h ^ has 
explicitly broken time-reversal symmetry, it is necessary 
to allow for the possibility of 3-spin couplings in the ef- 
fective model for the b spins. The effective Hamiltonian 
of the b spins £f e ff(5f,) is therefore of the form 

Hes(Sb) = —Jbb ^2 &bi&bj — Jbbb &bi(Tbj &bk ~ ^eff P" bi 

(23) 

which is an Ising model on the kagome lattice with three- 
spin interactions occuring within each 6-trimer. (A brief 
numerical study of a kagome Isin g m odel with three-spin 
interactions can be found in Ref. Il3l.) 

In the same manner in which we proceeded in 
Sec. 1111 Al we can integrate out the a spins to form an 
effective model for the b spins. 



Z = Tiexp(-/3H) 
= > Tr exp 

= £■ 

{S b } 



-PH{S b ,S a ,h) 



,-/3H el t(S b ) 



(24) 
(25) 

(26) 



In the presence of an applied field, time-reversal sym- 
metry is explicitly broken, and so there arc four distinct 
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terms in the trace. Z(S b i , S b 2, S b s) for each 6-trimer must 
be matched to the new effective model as follows (remem- 
bering to divide h c g by two because it is shared between 
two adjacent hexamers): 



m = l/3\ h eff- _2J bbb +4IJbb' 





= z a 


exp 


(3(3 J bb + 


Jbbb "+ 


3h c ff 

2 ' 


, (27) 




= z a 


exp 


P(-Jbb - 


~ Jbbb ' 


~2~ . 


, (28) 


Z(Ui,h) 


= z a 


exp 


P{-J bb H 


- Jbbb - 


hog, 
~ ~2~ . 


, (29) 


Z(iilh) 


= z a 


exp 


P(3Jbb - 


Jbbb - 


3/leffs 

2 ' 


• (30) 



From these equations, we find 

RT _i, z(m,h)z(ui,h) 
P Jbb ~ z ln z{in,h)z(ni,hy 

l^Z{m,h)Z(Ui,hf 



(3J t 



bbb 



(3hcS= 4 ln Z(Ul,h)Z(ltl),h> 
Z a = exp(f a ) = [Z(m,h)Z(Ui,h)] 1/8 

x [z(jn,h)z(m,h)f 8 . 



(31) 
(32) 

(33) 

(34) 
(35) 

The values of Z(S bl , S b2 , S b3 ) are given by Eqs. ([71)-([T0]). 
For h — > 0, we find that J bbb — > and h e g — > 0, and 
these equations reduce to our previous results. Eq. (|14|) ). 
In Fig. [3 we show how the external field h changes the 
effective coupling J bb as a function of temperature and 
Jaa- We see that a small field can change the sign of J bb 
when J^ n is weak. 



B. Ground state phase diagram of the 3-spin 
kagome Ising model 

We now have an exact mapping from the TKL in finite 
field to the kagome Ising model with the three-spin inter- 
action term J bbb . As we will see in the next section, the 
full ground state phase diagram of the original model will 
include phase transitions which only involve local a-spin 
configurations. However, we can develop some intuition 
about the long-range physics by studying the effective 
model for the b spins. Since this is now a classical model, 
the zero temperature phase diagram can be found in the 
usual manner by calculating the energy levels of each unit 
cell in order to determine the ground state: 

3 

^fe(TTT) = -3J&6 - Jbbb - ^^off, (36) 
£*(tti) = Jbb + Jbbb - -h eS , (37) 
Sfc(tU) = Jbb - Jbbb + \h eS , (38) 
E k (Ul) = -3J bb + J bbb + h cS . (39) 




m=l/3 


\ m=l 


t t 


\h e ff =- 2-bbb 


1 






\(0,0) h eff =0 


m=-l 


\ m=-l/3 


1 1 


\ 1 1 


1 



(b) J b , 



h eff=- 2 W+ 4IJ bbl X 




m=l 


m=l/3 t 

h e ff=2Jbbh 






(0,0) 


h eff= l2J bb l 

m=-l/3 


heff=-l24 1 \ 






m=-l | | \ 
1 




-2Jbbb- 4 l J bb' 



(c) J bb < 

FIG. 8: (Color online) Ground state phase diagram of kagome 
Ising model with three-spin interactions, (a) Ferromagnetic 
Jbb > 0. (b) Jbb = 0. (c) Antiferromagnetic Jbb < 0. m is the 
average magnetization per site. 



Representative ground state phase diagrams for this 
model are shown in Fig. [8] At large enough effective field 
strength \ h c g\, the system goes into a saturated ferromag- 
netic phase. However, large magnitude of the three-spin 
interaction \J bbb \ can frustrate this effect. In large parts 
of the phase diagram, we find that the ground state of a 
6-trimer is in the (til) state (or its time-reversed counter- 
part). When -Efc(HI) (-E'fe(tll)) is favored, then under- 
standing the character of the macroscopic ground state is 
reduced to the problem of enumerating the ways of tiling 
the kagome plane with one down (up) spin per kagome 
triangle. As was shown in Ref. [U, this is equivalent to 
placing dimers on the bonds of a honeycomb lattice. (See 
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Fig. 9 of Ref. 0-) In this phase, the correlation function 
(<7b(0)cr;,(r)) is equivalent to the dimer-dimer correlation 
function, which has been shown^ to be a power law. 1/r 2 . 

In the presence of the three-spin coupling, the topol- 
ogy of the phase diagram is dependent upon the sign of 
Jbb- When Jbb > 0, phase transitions between the two 
oppositely polarized dimer phases (i.e. from (ft!) to 
(Til)) are forbidden, and the dimer phases are separated 
by the saturated ferromagnetic phases. However, it is 
possible to go directly from one saturated ferromagnetic 
phase to its time-reversed counterpart. When Jbb < 0, 
then phase transitions directly from one dimer phase into 
its time-reversed counterpart are allowed, whereas phase 
transitions between the two saturated phases are not al- 
lowed. 



C. Full T = phase diagram of the XXZ-Ising 
model in finite magnetic field 

The full T = phase diagram can be obtained by 
adapting the method described in Sec. IV- A of Ref. 0- For 
every configuration of the three 6-spins in a hexamer, we 
have tabulated the energies of the eight 'internal' states 
of the hexamer in Table HI That is, we know the energies 
of the 64 states of the hexamer (including a- and 6-spins). 
The pattern of ground states within this 64-element ma- 
trix (for a given parameter set) allows us to determine the 
global ground state of the system.— Table Mil illustrates 
this procedure for a particular choice of parameters. 
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0.75 
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0.75 


0.75 


-2.3 


-2.3 


TTI 
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0.25 
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-1.8 


-1.2 
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-2.5 


0.98 


Til 


-0.25 


3.8 


-0.25 


-2.3 


-1.7 


1.2 


-0.9£ 


5 2.5 




3.2 


0.25 


0.25 


5.2 


-0.75 


1.2 


-1.8 


-1.8 



TABLE III: Energies of hexamer states for the parameters 
Jaa = J'aa = — 2, J at — — b h = 2 to two significant figures. 
The rows correspond to different (classical) configurations of 
the 6-spins (TTI represents the three configurations []], TITi 
and TTI), while the columns correspond to different (quan- 
tum) states of the a-spins. The ground state of the hexamer 
has energy —2.5 (indicated by the box), and it has the fa-spins 
in one of the three configurations |TTi TIT) or TTI- This im- 
plies that the ground state of the system is the honeycomb 
dimer phase (phase IV). 

In principle, we could, by comparing the analytic ex- 
pressions for the energy levels (Table [I]), derive analytic 
expressions for the phase boundaries. In this paper we 
will just present a phase diagram obtained by taking a 
grid of closely-spaced parameter points and finding the 
hexamer ground state(s) numerically for each parameter 
point. 

Fig. [§] shows the results of this procedure in the limit 
where the a-spins have isotropic Heiscnbcrg couplings 
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FIG. 9: (Color online) Phase diagrams of the TKL 
Heisenberg-Ising model for (a) the a-spins in the Heisenberg 
limit, J^a. = Jaa: an d (b) the a-spins in the Ising limit, 
Jaa = 0. In both panels, we use J ab = — 1 and T = 0. 



{Jaa = Jaa): as weu as m * ne previously studied case 
of Ising a-spins ( J^ a = 0) for comparison. Most of the 
phase boundaries in Fig. 9(a) appear to be straight lines, 
but four of them are slightly curved. The topology of 
Fig. 9(a) is more complicated than that of Fig. [51 because 



there may be distinct a-trimcr configurations underlying 
the same o-spin phase and because phase boundaries in 
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Fig. [8] may appear as extended regions of phases here. 
The "primed" phases are time-reversed versions of the 
"unprimed" phases. By examining the energy levels in 
detail (Table [TTT|) . we have identified the nature of each 
phase. 

Starting from the bottom of phase II and going through 
the phase diagram counterclockwise, the various phases 
are given as follows. Phase II is the 1(3/2) phase in a 
ferrimagnetic arrangement, so that the a-spins and 6- 
spins are perfectly aligned antiferromagnetically due to 
the coupling J* b , with the a-spins along the direction of 
the field h. With large enough h, the 6-spins become 
all up also, and there is a transition from phase II into 
phase I, which is a saturated ferromagnet. Increasing the 
strength of the a-trimer Hcisenberg coupling sufficiently 
one reaches phase X. In this phase, the 6-spins remain 
ferromagnetically aligned, but each a-trimer is now in a 
S = 1/2, S z = 1/2 state, which is doubly degenerate. 
The entropy per unit cell is therefore 2 In 2 = 1.3863.... 
Lowering the field takes us then into phase IV, which is a 
critical phase with power-law correlations that is equiv- 
alent to close-packed dimcrs on the honeycomb lattice 
(see Sec. HVBp . It is instructive to compare Fig. |9(b)| 
to Fig. 9(a) to see how the phase diagram evolves when 
the transverse coupling J^ a is introduced: phases I, II, 
and IV survive quantum fluctuations of the a-spins in- 
troduced by J^ a ^ 0, although the states of the a-trimer 
are somewhat modified. 

In the next section we give a detailed analysis of phase 
XL 
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TABLE IV: Energies of hexamer states, to two significant 
figures, for the parameters J,f a — J^ a = —2, J^b = —0.25, 
h — 2 (corresponding to phase XI, the "kagome-loop" phase). 
Boxes indicate ground states. 




FIG. 10: (Color online) An example of a configuration of b- 
spins in phase XI. Such configurations can be visualized as 
loops on the kagome lattice (see text). The statistical weight 
of the configuration is proportional to 2"', where I is the 
number of links in the loop. 



D. Phase XI: a "kagome loop gas" 

Phase XI is a new and unusual phase. A typical pat- 
tern of hexamer energies for phase XI is shown in Ta- 
ble IIVI This means that in phase XI, each 6-trimer is 
allowed to contain either down spins or 2 down spins, 
but the latter case is less likely by a factor of 2. The con- 
straint that each 6-trimer must contain an even number of 
down-spins means that the down-spins in the lattice must 
form nonintersecting closed loops - i.e., allowed configu- 
rations correspond to Eulerian subgraphs of the kagome 
lattice. (An Eulerian subgraph is one in which the degree 
of every vertex is even.) However, not every Eulerian sub- 
graph corresponds to an allowed configuration of phase 
XI, because we are not allowed to put three down-spins 
on the same 6-trimer. Figure \W\ shows an example of 
an allowed configuration. In moving through the "phase 
space" of the ground state manifold, loop number is not 
conserved, loops are non- intersecting (i.e. they are re- 
pulsive with a hard-core interaction), and loop size is 
not conserved (i.e. the loop tension is zero). This cor- 
responds to a non-intersecting loop gas on the kagome 
lattice. 

The entropy in phase XI is non-trivial. First, we 
present an estimate based on counting loops on the 
kagome lattice. This is a series expansion using a; = 



= i as a small parameter. Consider a system 
with V unit cells. The most probable configuration of 
the 6-spins is the one with all spins up (Fig. 11(a) I. The 



next most probable state has a hexagon of down-spins 
(Fig. 1 1 (b) I . This state has the same energy (the ground 



state energy) , and comes with a factor of V because the 
hexagon can be placed anywhere in the lattice, but its 
statistical weight is smaller by a factor of x 6 , because 
there are 6 triangles that are (til) instead of (TTT) - Con- 
tinuing in this fashion (Fig. |ll(c)|ll(d)|ll(e)[ ), we obtain 
a series for the total number of states: 



W 



-,2V 



1 + Vx b + 3Vx L 



2V + V{V - 7y )^ 



so the entropy per unit cell is 
In IF 



lim 



V 



2 In 2 



3a; 



10 



3.r 



12 



(40) 



(41) 



i, including progressively more terms in 



When x 

the series gives the approximations 1.38629; 1.40192; 
1.40485; and 1.40448, suggesting that the entropy per 
unit cell is s = 1.405(1). 



11 




FIG. 11: (Color online) Diagrams representing terms of up to 
order x 12 in the series for the entropy of phase XI (when h is 
small and positive) . Gray lines represent a background of up- 
spins on the kagome lattice. Bold (red) lines represent loops 
of down-spins. Shaded triangles represent triangles of fa-spins 
in an |J,1 configuration; each shaded triangle corresponds to 
a factor of x = \ in the statistical weight. 



The mean magnetization of the 6-spins can be calcu- 
lated by a similar series expansion: 



To;, = 3 — 2 (6a 



30a; 10 - 18x 12 



(42) 



When x = \ this series gives m& as 0.12(1). It seems 
likely that the true result is an irrational number. This is 
in contrast with the honeycomb dimer phase (phase IV), 
which has rational magnetizations on both sublattices, 
m a = m b = 1/3. 

We now present an exact calculation of the entropy of 
phase XL First, note that the kagome lattice is related to 
a honeycomb lattice in the following way. Every triangle 
of the kagome lattice maps to a site of the honeycomb lat- 
tice, and every site (spin) of the kagome lattice maps to a 
bond on the honeycomb lattice. Recall that in phase XI, 
every triangle on the kagome lattice contains cither or 
2 down spins (call these "active sites"). Therefore, every 
allowed spin configuration on the kagome lattice corre- 
sponds to an "active-bond" configuration on the honey- 
comb lattice, in which every site of the corresponding 
honeycomb lattice must have either or 2 active bonds 
adjacent to it. (Sec Fig. [T2l) 

The higher-order terms in the series for the degeneracy 
of phase XI (Eq. ([4T)]) ) can be seen to be identical to the 
high-temperature series for the partition function of the 
Ising model on the honeycomb lattice: 



phoney 
^Ising 



(cosh/3J Is ; 



honey \3V2^ 
sing / 



1 + vt 6 + wt 
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FIG. 12: (Color online) Illustration of the mapping from a 
configuration of loops on the kagome lattice (bold red lines) 
to loops on a honeycomb lattice (dashed blue lines). 



when t = tanh (3 J —1/2 and V is the number of honey- 
comb unit cells. The entropy of the kagome loop phase 
is similarly related to the free energy of the honeycomb 
Ising model. Using the well-known exact integral expres- 
sion for the latte r 10 ! 11 , we find that the entropy of phase 
XI (per unit cell) is 



TKL _ /•honey 

^ phase XI Rising 

= 1.4053. 



31ncosh/3Jj h ™ g cy 



(44) 



(43) 



which agrees with the series result presented earlier. 

It is interesting that the kagome loop gas at T = 
and h ^ can be exactly mapped to a honeycomb Ising 
model at T > and h = 0. 



V. CONCLUSIONS 

In conclusion, we have studied an XXZ-Ising model on 
the triangular kagome lattice, in which spins residing on 
small triangles (a-trimers) are fully quantum mechani- 
cal, and spins residing on large triangles (6-trimers) are 
classical Ising spins. For all fields and temperatures, we 
have shown that there is an exact mapping to the kagome 
Ising model. 

In the absence of applied field, the full Ising limit 
(i.e. J^a = 0) of the model has a phase transition 
from an ordered phase (Phase I) for J^ a above a crit- 
ical value Jad mt = —\Jab\-> to a disordered phase with 
f° r Jaa < — \Jab\ wm ch nas residual entropy per unit 
cell ln72 = 4.2767... (Phase V). The introduction of an- 
tifcrromagnctic quantum fluctuations of the a-spins (i.e. 
J'aa < 0) nas little effect on the ordered phase, but it has 
dramatic effect on the disordered phase, partially lifting 
the degeneracy of the ground state and reducing the en- 
tropy by almost a factor of two to 2.5258... per unit cell. 
For strong ferromagnetic quantum fluctuations of the a- 
spins (i.e. J^ a > 0), a new phase (phase 7(1/2)) appears 
in which the a-trimers have S = 3/2; S z = 1/2. There 
is a zero-temperature phase boundary on the XY side of 
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the Heiscnbcrg line, for J* a = J z a + l^abl- separating this 
phase from the one with fully polarized a-trimers, phase 
7(3/2), with S = 3/2; S z = 1/2 within each a-trimcr. 

All of the ground state phase transitions which occur 
in zero field survive at finite temperature, except for the 
ground state transition from an ordered phase with S = 
3/2; S z = 3/2 (i.e. Phase 1(3/2)) to an ordered phase 
with S = 3/2; S z = 1/2 (i.e. Phase 7(1/2)). Since this 
transition is entirely about local physics, it becomes a 
crossover at finite temperature. However, its effects are 
quite evident at finite temperature, as the 7(3/2) phase 
is more robust, and survives to much higher temperature 
than does the 7(1/2) phase. 

In the presence of applied field, we find that the XXZ- 
Ising model on the TKL maps to the kagomc Ising model 
in field with 3-spin interactions. We find exact solutions 
for the ground state of this model. In the absence of 
quantum fluctuations, but for strongly frustrated inter- 
actions, we have previously shown that the Ising TKL has 
a critical spin liquid ground state (with power-law cor- 
relations) for any weak applied field (Phase IV), which 
can be seen through a mapping to hard core dimcr cov- 
erings of the honeycomb lattice. We find that this phase 
survives the addition of quantum fluctuations on the a- 
trimers; however, it is then necessary to apply a finite 



magnetic field h to access this phase (rather than an in- 
finitesimal field, as in the case of the TKL Ising AF). 

For antifcrromagnctic interactions, we have shown that 
small (and even infinitesimal) applied fields lead to a 
"kagome loop gas" ground state (phase XI). We have 
mapped this to the ferromagnetic honeycomb Ising model 
at zero field and finite temperature, and hence we have 
shown that phase XI has an entropy of 1.4053... per unit 
cell. It is likely that the &-sublatticc magnetization is ir- 
rational, in contrast to all the other phases studied in the 
paper. 
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In general, knowledge of the finite-temperature thermody- 
namics of a subsystem does not necessarily allow one to 
deduce the thermodynamic behavior of the entire system. 
It is fortunate that the TKL is simple enough such that 
knowledge of the ground state of a single hexamer allows 
one to infer the ground state of the entire system. 



